Luttinger parameters and momentum distribution function for the 
half-filled spinless fermion Holstein model: A DMRG approach 
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Abstract. - We reexamine the nature of the metallic phase of the one-dimensional half- filled 
Holstein model of spinless fermions. To this end we determine the Tomonaga-Luttinger-liquid 
correlation parameter K p by large-scale density-matrix renormalisation-group (DMRG) calcula- 
tions, exploiting (i) the leading-order scaling relations between the ground-state energy and the 
single-particle excitation gap and (ii) the static charge structure factor in the long-wavelength 
limit. While both approaches give almost identical results for intermediate-to-large phonon fre- 
quencies, we find contrasting behaviour in the adiabatic regime: (i) K p > 1 (attractive) versus (ii) 
K p < 1 (repulsive) . The latter result for the correlation exponent is corroborated by data obtained 
for the momentum distribution function n(k), which puts the existence of an attractive metallic 
state in the spinless fermion Holstein model into question. We conclude that the scaling relation 
must be modified in the presence of electron-phonon interactions with noticeable retardation. 



During the last three decades we have seen a constant 
growth of experimental realizations of one-dimensional 
(ID) materials. Nowadays the progress in nanotechnol- 
ogy allows to manufacture isolated carbon nanotubes or 
quantum wires [1] . But there are also bulk materials with 
quasi-lD structures inside. Famous examples are con- 
jugated polymers, charge transfer salts, halogen-bridged 
transition metal complexes, ferroelectric perovskites, spin 
Peierls compunds, molecular metals or organic supercon- 
ductors [2]. The apparent diversity of physical proper- 
ties observed for different material classes has its seeds in 
the strong competition between the itinerancy of the elec- 
tronic charge carriers on the one hand and the electron- 
electron and electron-lattice interactions on the other 
hand. The latter tend to establish insulating spin-density- 
wave or charge-density-wavc (CDW) ground states, re- 
spectively, at least for commensurate band fillings [3-8]. 
Interactions have drastic effects in ID systems compared 
to higher dimensions. Most notably one observes a "col- 
lectivisation" of any excitation. As a consequence, for 
fcrmionic systems, the usual Fermi liquid description 
breaks down [9] . Luttinger liquid theory provides an ade- 
quate compensation [10]. It tells us that all ground-state, 
spectral and thermodynamic properties of a Luttinger liq- 
uid are basically controlled by a few (non-universal) pa- 



rameters. This result can be used in the following way. 
Starting from a specific microscopic model, one can try to 
compute certain (thermodynamic) quantities exactly, e.g. 
for finite systems by elaborate numerical techniques, and 
afterwards extract the Luttinger liquid parameters, e.g. 
the charge correlation exponent K p and charge velocity u p , 
out of them. Advantageously these parameters, describing 
the overall low-energy physics of our system, are much less 
sensitive to finite-size effects than the correlation functions 
themselves. Of course, the concept of a Luttinger liquid 
has to be taken as a starting point to study more complex 
situations, comprising e.g. the lattice degrees of freedom 
or disorder effects. Then the Luttinger parameters be- 
come effective parameters, which characterise very basic 
properties of the system, such as an attractive {K p > 1) or 
repulsive (K p > 1) interaction between the particles [9,11]. 

In this respect, focusing on the coupling of charge carri- 
ers to the vibrations of a deformable lattice, the so-called 
Holstein model of spinless fermions (HMSF) [12], 
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is particularly rewarding to study. It accounts for a tight- 
binding electron band (oat), a local electron-phonon (EP) 
interaction (oc g), and the energy of the phonon subsys- 
tem in harmonic approximation. In eq. (1), cj (c.) creates 
(annihilates) an electron at Wannicr site j of a ID lattice 
with N sites, b] (bj) are the corresponding bosonic op- 
erators for a dispersionlcss optical phonon, and loq is the 
frequency of such an internal (e.g. molecular) vibration. 

Despite its seemingly simplicity, the ID HMSF is not ex- 
actly solvable. It is generally accepted, however, that the 
model exhibits a quantum phase transition from a metal 
to a CDW insulator at half-filling, when the EP coupling g 
increases at fixed ujq > 0. 1 The CDW phase above g c (u)o) 
is connected to a (Peierls) distortion of the lattice, and 
can be classified as traditional band insulator and pola- 
ronic supcrlattice, respectively, in the adiabatic (loq <C 1) 
and anti-adiabatic (u>q 3> 1) regimes [13,14]. A wide range 
of analytical and numerical methods have been applied to 
map out the phase diagram of the HMSF in the whole 
g — u>o plane [13, 15-19], with significant differences in the 
region of small-to- intermediate phonon frequencies. The 
results agree in the anti-adiabatic strong EP coupling limit 
(ujq — > oo, g > 1), where the HMSF possesses XXZ-modcl 
physics. There a Kostcrlitz-Thouless type transition [20] 
occurs at the spin isotropy point, with K p reaching 1/2 
from above at the transition point [15, 18]. 

In the first instance, however, the correlation expo- 
nent K p can be used to characterise the metallic phase 
itself. According to Haldane's conjecture [21], a ID gap- 
less (metallic) system of interacting fcrmions should be- 
long to the Tomonaga-Luttingcr liquid (TLL) universality 
class [10,22]. For a TLL of spinless fermions, the ground- 
state energy E(N) and the one-particle charge excitation 
gap A Cl of a finite system with TV sites scale to leading 



order as [11,23]: 



E(N) « P 1 

= e(oo) — -, 

N V ' 3 2 /V 2 ' 

A C1 = E ± (N) - E(N) = — — . 

V ' K ' 2 K p N 
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Here, e(oo) denotes the energy density of the infinite sys- 
tem with N/2 electrons, E ± (N) are the ground-state en- 
ergies with ±1 fermion away from half filling, and u p is 
the renormalised charge velocity. 2 Using these equations 
from field theory, in the past, K p and u p were determined 
for the HMSF on finite clusters by various exact numeri- 
cal techniques [13, 17, 18,26]. Interestingly the TLL phase 
seems to split into two different regions: For large phonon 
frequencies the effective fcrmion-fermion interaction is re- 
pulsive (K p < 1), while it is attractive (K p > 1) for small 
frequencies [13,26]. In the former (anti-adiabatic) regime 



We consider the half-filled band case hereafter, i.e., 
i (ct c . ) = ^, and take t = 1 as energy unit. 

2 The TLL scaling relations (2) and (3) were also derived for spin- 
ful systems [24], and e.g. used in order to compute the central charge 
in the framework of the t—J model [25] . 
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Fig. 1: (Colour online) Finite-size scaling of the ground-state 
energy E(N) and the one-particle charge excitation gap A C1 
(inset) in the spinless Holstein model at half filling. Results 
are obtained by DMRG for c^o =0.1 and g — 2.5. The linear 
equations give the coefficients of a straight-line fit to the scaling 
relations (2) and (3). 



the kinetic energy (oc u p ) is strongly reduced and the 
charge carriers behave like (small) polarons [27,28]. By 
contrast the mass renormalisation is rather weak in the 
adiabatic regime [26]. The size of the phonon frequency 
also significantly affects the electron and phonon spectral 
functions [14,19,29]. 

The existence of an attractive TLL (K p > 1) in the 
HMSF is by no means obvious however. Although retar- 
dation effects might lead to an effective attraction between 
electrons at small oJo/t (i.e. , a second electron may take the 
advantage of the lattice distortion left by the first one), it 
has been pointed out that such an interaction is ineffective 
in the case of spinless fermions for small EP couplings be- 
cause of the Pauli exclusion principle [15]. 3 Furthermore, 
if K p would increase with increasing EP coupling at small 
luq, as indicated by different numerical studies [13,26] ex- 
ploiting equations (2) and (3), how could we detect the 
phase transition from K p — > 1/2 in the adiabatic regime? 
Of course, equations (2) and (3) are leading-order expres- 
sions only, and nonlinear correction terms have to been 
taken into account in order to obtain accurate data for 
g c [18]. This particularly applies to the adiabatic region. 
According to table III in Ref. [18] the difference between 
the g c determined with and without nonlinear correction 
terms adds up to more than 3% for loq = 0.1, whereas it 
is only 0.4% for ujq = 10. Actually the charge velocity 
Up depends strongly on the system's size as shown in fig- 
ure 1 obtained by a density matrix renormalisation group 
(DMRG) calculation. Extrapolating the ground-state en- 
ergies E(N) for N = 8, 12 and 16, the charge velocity u p 
can be estimated as u p /2 ~ 0.977, so that K p ~ 1.21 from 
the finite-size scaling of A Cl , while taking the ground-state 
energies for N = 32, 48 and 64, the extracted value of K p 
reduces to 1.06 (u p /2 ~ 0.858). Moreover we are faced 



'Note this argument does not hold for the spinful Holstein model. 
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with the difficulty that the single particle excitation gap 
seems to scale to zero (see inset of figure f for uj = 0.1, 
g — 2.5), i.e. A c i gives no signal for a pairing instability. 
Because of this situation it is highly desirable to find a 
reliable and numerical efficient method for calculating the 
correlation exponent K p with high precision in the whole 
TLL regime. 

Recently Ejima et al. [30] have shown that K p can be 
determined for fermionic models accurately in an alterna- 
tive way: By a DMRG calculation of the charge structure 
factor for systems with open boundary conditions. The 
approach was extended to coupled fermion-boson systems 
and has been used, e.g., to analyse the metal-insulator 
transition points in a model with boson affected transport, 
for both small and large boson frequencies [31]. 

In this work we adapt this calculation scheme to the ID 
Holstein model of spinless fermions (1) and reexamine, in 
particular, the possible existence of a metallic phase with 
attractive interaction. To this end we compute, in a first 
step, the static charge structure factor 
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and extract, in a second step, the TLL correlation expo- 
nent K p , being proportional to the slope of S c {q) in the 
long-wavelength limit [7,30,32]: 

Sc (q) 



K D = 7r lim 
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Moreover we calculate the momentum distribution func- 
tion for the HMSF and, having accurate data for K p at 
hand, analyse the results within a TLL description, also 
in relation to the corresponding results for the half-filled 
spinless t-V model. When treating coupled ID fermion- 
boson systems by DMRG we employ the pseudo-site ap- 
proach [33] which maps a bosonic site, containing 2 nb 
states, exactly to rib pseudo-sites. For the numerics pre- 
sented below we have taken into account up to rib = 5 
pseudo-sites, so that the ra^-th local boson density is al- 
ways smaller than 10~ 8 . In addition we kept m = 1200 
density-matrix eigenstates. Then the discarded weight 
was always smaller than 1.0 x 10 -10 . 

Figure 2 (a) presents K p obtained from eq. (5) for var- 
ious phonon frequencies. Note that data points in (a) 
represent K p - values extrapolated to the infinite system at 
fixed (g, u)q) [cf., e.g., panel (c)]. For intermediate-to-large 
phonon frequencies we find K p (g) < 1 for all g, but an 
appreciable reduction of K p takes place above 3 = 1 only. 
The strong decrease of K p and u p (not shown) for g > 1 
is closely connected to polaron formation, which appears 
at about g ~ 1 in the non-to-anti-adiabatic regime [28]. 
There the TLL typifies a (repulsive) polaronic metal [14]. 
We emphasise that in this frequency region the values of 
A' p , computed from eq. (5) via the static charge struc- 
ture factor, reasonably agree with those determined by 
the scaling relations (2), (3) [see panel (a), filled and open 
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Fig. 2: (Colour online) Panel (a): TLL parameter K p in the 
spinless Holstein model at half filling. Closed symbols are 
obtained via S c (q) from eq. (5) for ljq = 0.1 (triangles), 1 
(squares), and 10 (circles). K p obtained from the scaling 
relations (2) and (3) are included as open symbols. Lines 
are guides to the eye. Panel (b): Ground-state phase dia- 
gram of the ID half-filled spinless Holstein model according to 
Refs. [13,14,18]. Stars denote the phase transition points ob- 
tained from K p = 1/2 in (a). Panel (c): K p as a function of 
the inverse system size at various EP couplings g = 1.5, 2.5, 
and 2.81 for uoa = 0.1 (adiabatic regime). Lines are polynomial 
fits. 



symbols for luq = 10]. Furthermore, our values for the 
critical coupling, g c , confirm previous results (although a 
Kostcrlitz-Thouless transition is difficult to detect because 
the gap opens exponentially slow), as can be seen by in- 
serting the points where K p (g c ) = 1/2 (stars) into the 
existing phase diagram of the HMSF [13, 14, 18] [cf. panel 
(b)]. 

Let us now look whether the situation changes when 
the phonon frequency becomes smaller, i.e., when we en- 
ter the adiabatic regime. Figure 2 (c) shows the scaling 
of K p at ojq = 0.1, based on the relation (5), for up to 
TV = 256 sites, using open boundary conditions. The lines 
are second-order polynomial fits. Surprisingly, we find 
that K p scales to values smaller than unity for any EP 
coupling (filled symbols). This holds for other adiabatic 
phonon frequencies ujq < 1 as well. Taking this result seri- 
ously, we arrive at the conclusion that the HMSF does not 
exhibit a metallic TLL phase with attractive interaction, 
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which is in strong contradiction to the reasonings based 
on the leading-order energy scaling laws (2) and (3) [see 
open symbols in figure 2 (a)]. We would like to point out, 
however, that our S , c (q)-bascd approach gives apparently 
the correct value of the critical coupling for the TLL-CDW 
metal- insulator transition in the adiabatic HMSF. In pre- 
vious work, g c was estimated as g c (too = 0.1) ~ 2.8 [18], 
which is in accordance with our DMRG- results for K p {N) 
at g = 2.81 that clearly extrapolate to K p = 1/2 in the 
thermodynamic limit. This means that the TLL-CDW 
transition at small uiq = 0.1 seems to be of Kosterlitz- 
Thouless type as well. 

To substantiate these findings, we investigate another 
quantity of interest, the so-called momentum distribution 
function, 







(G) 



Basically n(k) is the Fourier transform of the equal time 
Green's function [9] and therefore gives the occupation of 
fcrmionic states carrying momentum k. For free fermions, 
at T = 0, all states up to the Fermi energy, Ep, are occu- 
pied, so that n(k) has a discontinuity (Z = 1) at the cor- 
responding Fermi momentum kp. A For a ID TLL, instead 
of the (Fermi liquid archetypical) jump of n(k) at kp, one 
finds an essential power law singularity, corresponding to 
a vanishing quasiparticle weight Z = 0, 



i{k) 



- C\k - k F \ a sgn(k - k F ) , 
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where nt F = 1/2 for the half-filled band case. For spinlcss 
fermions, again the critical exponent a is given by the TLL 
parameter K p : 



1 
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The relation (7) with (8) was first derived in [10,34], and 
afterwards many analytical [35] and numerical [36] cal- 
culations were performed in order to determine the mo- 
mentum distribution in the weak and strong coupling 
regimes. By means of DMRG, n(k) can be computed 
directly from the Fourier transformed (c^c ; ) correlator, 

n ( k ) = 7?J2f,i=i C0S ( k (j ~ 0) (c} c i) ' where k = W m 
with m = 0, ... , N/2. In the following, we calculate n(k) 

for a linear chain with periodic boundary conditions, and 

TV = 66 sites. 

Before we discuss n{k) for the HMSF, let us consider a 
somewhat simpler, purely fcrmionic model, however, the 
spinlcss t-V model, 



H = -t£ 



(9) 



where V is the nearest neighbour Coulomb interaction. 
This is of avail because the t-V model can also be mapped 



4 In an interacting Fermi liquid system there is still a discontinuity, 
but Z < 1. 



1c 






0.8 




(a) 


0.6 

c; 

0.4 


K = 0.75 

p 


- 

V= 1 


0.2 


- C= 0.515 


(TLL) 




1c 




< ^Oe>r=H°Aa .-. -i 

gj^J ivu UUUUOO O Ql 






r ouQ aasoo6o„' 




0.8 




(c); 


0.6 

c 

0.4 




- 


K = 0.556 


V=-\.9 


0.2 


p 

■ C=0.55 


K (TLL) 






, °?°5qoBorx». 


) 0.2 0.4 


0.6 0.8 



K = 0.649 

p 

C=0.53 



(b) 



V= 1.5 

(TLL) 



oooooooo 



(d) 



V= 1000 

(CDW) 



0.2 



kin 



(CDW) 

0^4 0^6 

kin 



O °°oooooo 



0.8 



Fig. 3: (Colour online) Momentum distribution function n(k) 
for the half-filled spinless t-V model. Lines are fit to eq. (7) 
with K p taken from eq. (10). 



onto the exactly solvable XXZ-Heisenberg model (i.e., 
it should exhibit the same asymptotic behaviour as the 
strong-coupling anti-adiabatic HMSF). For the t-V model 
the analytical form of K p in the thermodynamic limit is 
known, 



X = 

p 2arccos[-y/(2t)] ' 



(10) 



Hence, the results obtained for n(k) by DMRG can be fit- 
ted by the relation (7), with K p taken from eq. (10) [37]. 
Clearly, since eq. (7) is a weak-coupling result, the DMRG 
data for n{k) are fitted almost perfectly for small V. 
This is demonstrated by figures 3 (a) and (b). Figure 3 
(c) shows that the agreement becomes worse for larger 
Coulomb interaction (V = 1.9). In the insulating phase 
(V > 2), the power-law singularity does not exist any- 
more [c.f. the smooth curves in figure 3 (d)]. As V — > oo, 
the system becomes a "perfect" CDW, and consequently 
n(k) = 1/2 for all momenta k. De facto this situation 
is realized for V = 1000 already where, according to fig- 
ure 3 (d), n(k) is almost constant. 

Turning now to the HMSF, we discuss at first the case of 
large phonon frequencies. Figure 4 gives the n(k) DMRG 
data obtained for ujq = 10 (symbols). Obviously, the mo- 
mentum distribution is a monotonously decreasing func- 
tion as k changes from the centre (k = 0) to the boundary 
of the Brillouin zone (h = 7r), with a power-law singular- 
ity at k = k-p in the metallic phase [panels (a) to (c)]. 
Quite different from the t-V model, however, the momen- 
tum distribution becomes renormaliscd for all momenta 
k, as soon as the EP coupling is switched on, where n(k) 
decreases (increases) almost uniformly for < k < ir/2 
(tt/2 < k < 7r). Although there is no jump in n(k) at kp, 
as for an ordinary Fermi liquid, for finite TLL systems 
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Fig. 4: (Colour online) Momentum distribution function n(k) 
in the anti-adiabatic regime (cjq — 10) of the half-filled spinless 
Holstein model . Lines are a fit to eq. (7) with K p calculated by 
the DMRG. Insets give n(k) with magnified axis of ordinate. 



Fig. 5: (Colour online) Momentum distribution function n(k) 
in the adiabatic regime (cjq — 0.1) of the half-filled spinless 
Holstein model. Lines are a fit to eq. (7) with K p calculated 
by the DMRG. 



the difference A = nikp — 6) — n(k-p + 5) is finite (with 
S = n/N = in our case), and rapidly decreases with 
increasing EP interaction g. As can be seen from the solid 
lines in figure 4 (a)-(c), the momentum distribution can 
be surprisingly well fitted to the weak-coupling result (7), 
just by adjusting the constant C. Thereby we take the 
ifp-values extracted from equation (5). Of course, around 
k — kp the agreement becomes worse as g increases, but 
we observe a power-law singularity even close to the CDW 
transition point. Approaching the insulating CDW state 
this singularity vanishes, and A — > as g — > g c [cf. the in- 
sets in panels (c) and (d)]. In the CDW phase, n(k) ~ 1/2 
for all k [sec panel (d)]. In the anti-adiabatic regime, the 
CDW state basically constitutes a polaronic superlattice, 
i.e the electrons are heavily dressed by phonons and, in 
addition, ordered in a A-B-structure. Since the polarons 
are self-trapped, the system tends to be a perfect CDW, 
as in the limit V — > oo of the t-V model. 

Finally, we investigate the behaviour of n(k) in the adi- 
abatic regime of the HMSF (see figure 5 for oj = 0.1). 
In this case, n(k) is well approximated by equation (7) 
with C = 0.5 for all g < g c [see panels (a) and (b)]. This 
means the weak-coupling result K p < 1, A < 1, holds 
in (almost) the entire metallic region, where the system 
can be considered as to be composed of nearly free elec- 
trons. The momentum distribution starts to deviate from 
equation (7) just in the neighbourhood of the metal in- 
sulator transition point g ~ g c , but even there n(k) dif- 
fers near k ~ kp — n/2 only. Note that A (C) is still 
very large in the transition region. Of course, very close 
to the critical point, where a strong renormalisation of 
K p takes place (indicating the formation of a TLL with 



strong repulsive interactions), the fit of our DMRG data 
to the weak-coupling relation (7) fails. In the insulating 
state, n(k) is given by a smooth curve (without power- 
law singularity), which - in contrast to the anti-adiabatic 
case - exhibits a significant curvature because the EP cou- 
pling used in figure 5 (d) is small if compared to the half 
electronic bandwidth 2t. 5 Therefore system now typifies 
rather a Peierls band insulator than a polaronic superlat- 
tice. 

To summarise, we investigated the properties of the 
metallic phase and the metal insulator transition in the 
spinless fcrmion Holstein model by means of a boson 
pseudo-site DMRG technique supplemented by a careful 
finite-size scaling analysis. In particular we determined 
the Tomonaga-Luttinger correlation exponent K p from 
the long-wavelength limit of the static charge structure 
factor. This approach yields reliable data for K p in the 
whole range of clcctron-phonon interaction strengths g and 
phonon frequencies u>q. We compare our results with new 
and previous data extracted in an alternative way from 
leading-order scaling relations for the ground-state energy 
and single-particle excitation gap. In striking contrast to 
the latter data we find K p < 1 for all phonon frequen- 
cies, i.e., the metallic state of the HMSF represents a re- 
pulsive Tomonaga-Luttinger liquid, even in the adiabatic 
regime. Therefore we conclude that in one dimension we 
have to include the spin degrees of freedom [7, 8, 15] in 
order to obtain e.g. a phase with attractive interactions 
{K p > 1), or even dominant superconducting correlations. 

5 Recall that in the adiabatic regime A = g 2 uJo/2t is the appropri- 
ate dimcnsionless EP interaction parameter in order to discriminate 
weak (A 1) and strong coupling (A 3> 1) situations. 
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Furthermore, since the metal insulator phase boundary in 
the g-ojQ 1 plane obtained from the K p (g,uJo) = 1/2 line 
is in excellent agreement with previous results [13,14,18], 
we suppose that the TLL-CDW transition in the HMSF 
is always of Kosterlitz-Thouless type. Comparing the be- 
haviour of the momentum distribution function with the 
weak-coupling TLL result reveals, however, significant dif- 
ferences regarding the nature of the metallic and insulat- 
ing phases in the adiabatic and anti-adiabatic regimes of 
the HMSF. Whereas the metallic state is a weakly renor- 
maliscd TLL and the CDW phase typifies a Peierls band 
insulator at small phonon frequencies, a polaronic metal 
and a polaronic superlattice are formed at large phonon 
frequencies. This is in accord with the electron and 
phonon spectral properties detected in Refs. [14,29]. In 
the strong-coupling anti-adiabatic regime the momentum 
distribution function indicates perfect CDW behaviour as 
in the V — > oo limit of the t-V model. 
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